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Abstract 

Discrete Dirac type self-adjoint system is equivalent to the block 
Szego recurrence. Representation of the fundamental solution is ob- 
tained, inverse problems on the interval and semiaxis are solved. A 
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block Toeplitz matrices are treated. 
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1 Introduction 



Continuous self-adjoint Dirac type system 



dY_ 

dx 



(x,z) = i(zj+jV(x))Y(x,z), j 







V 



v 
v* 



is a classical object of analysis with various applications (in mathematical 
physics and nonlinear integrable equations, in particular). Here I p is the pxp 
identity matrix and v is a p x p matrix function. In this paper, we treat a 
discrete self-adjoint Dirac type system: 



W k+1 {X)-W k (X) 



-jjC k W k (X) (k>0), 



;i.2) 



i 



where C& are m x m matrices, m = 2p, 



Ck — Cfc, CkjCk — j, 



(1.3) 



To see that (11.21) is a discrete analog of system (11.11) . notice that (11.11) 
is equivalent to a subclass of canonical systems W x = izjH(x)W (see [331 
137] and references therein). One can follow also the arguments from |24J, 
where the skew self-adjoint discrete Dirac type system have been studied 
and explicit solutions of the isotropic Heisenberg magnet model have been 
obtained. As suggested in |24J introduce matrix functions U and W by the 
relations 



W(x,z) = U(x)Y(x,z), —(x) = -iU(x)jV(x), 17(0) = I m . (1.4) 



Since V is self-adjoint, we get from ( 11.41) that U is j-unitary, i.e., UjU* = j. 
Now (II. ip and the first relation in (11.41) yield 



where H = jUU*j = H*, HjH = j. Compare system (11.21) . where matrices 
Ck satisfy (11.31) . and system (II. 5p to see that (11.21) is an immediate discrete 
analog of (11.11) . 

When p — 1 and Ck > 0, system (II. 2p is equivalent to the well-known 
self-adjoint Szego recurrence, which plays an important role in the orthogonal 
polynomials theory and is also an auxiliary system for the Ablowitz-Ladik 
hierarchy (see, for instance, P33 [2U1 EB] and various references therein). The 
equivalence of system (ll.2p . where Ck > and CkjCk = j, to the block 
(matrix- valued) Szego recurrence is given in Proposition 12.11 

We consider representation of the fundamental solution of system (11.21) 
and solve direct and inverse problems directly in terms of the Weyl functions. 
Both explicit and general solutions are obtained. First, we obtain explicit so- 
lutions of the direct and inverse problems for system (II. 2p for the case of the 
so called pseudo-exponential potentials Ck (the case of the rational Weyl func- 
tions). Our case includes as a subcase the rapidly decaying strictly pseudo- 
exponential potentials. Recall that discrete and continuous systems with the 
potentials, which belong to the subclass of the strictly pseudo-exponential 



dW 



dU_ 

dx 



U- X W + iU(zj + ]V)U~ l W = izjHW, 



(1.5) 



dx 
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potentials, have been actively studied in PQ-JH], [E]-[in]- In particular, direct 
and inverse problems for Szego recurrence on the semiaxis with the scalar 
(p = 1) strictly pseudo-exponential potentials have been treated in [51 [6]. 
Direct and inverse problems for the pseudo-exponential potentials (contin- 
uous case) have been studied in a series of Gohberg-Kaashoek-Sakhnovich 
papers (2TJ [22] (see references therein and see also [T7| for the case of the 
generalized pseudo-exponential potentials). The case of the discrete skew- 
self-adjoint Dirac system have been studied in [21]. Notice that similar to 
[2H [31] (see also [5], [6], [HI [2TJ [23] ) we start our explicit constructions with the 
explicit formula for the fundamental solution. 

For a more general (non-rational) situation of the Weyl functions <fi(z) = 



such that J2T=o ll^fcll < 00 ^ ne direct problem for a block (matrix- 

k=0 

valued) Szego recurrence on the semiaxis (including non-self-adjoint case and 
under some additional conditions) is treated in a recent important Alpay- 
Gohberg paper [Tj. In Sections and M we solve direct and inverse prob- 
lems for the general type potentials Cj. > (and thus for the general type 
self-adjoint block Szego recurrence) on the interval and semiaxis. Borg- 
Marchenko type uniqueness result for system (jl.2p is obtained too. Con- 
nections with the well known Toeplitz matrices appear. For the interesting 
discussions on the connections between Toepliz matrices, Szego recurrencies 
and orthogonal polynomials see also J3J [T3] and references therein. Inter- 
esting spectral theoretical results on the discrete canonical systems, where 
CkjCk = 0, one can find in [271 EE]- A complete Weyl theory for the Jacobi 
matrices and various useful references are contained in [39] . 

2 Preliminaries 

An important discrete analog of Dirac type system takes the form 



CO 




Xk+i(z) — k R k 



L 



Xk(z), Rk = Rl, RkjRl = 3 



(2.1) 



v J 
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(9 k E C, 9k 7^ 0). System (12. ip can be presented in the form f 1 1.2 1) after 
transformation 

Xk(z), C k = (Uir 1 RlU k -\ (2.2) 
where U := I m , 

U k :=(ijR )(ijR 1 )x...x(ijR k _ 1 ) (k > 0), 2 = (2.3) 

1 — 1\ 

A particular scalar case (p = 1) of system (12. ip is a well known Szego recur- 
rence, where 

>0, |p fc |<l, 9 k = ^l-IPfel 2 . (2.4) 

Further we assume that some expression for 9k via is fixed for p > 1 too. 
For p = 1 representations (12.41) of follow from the relations Rk = R* k > 
and RkjRk = j- Coefficients pk are called Schur (or sometimes Verblunsky) 
coefficients (see, for instance, [T2l [38] and various references therein). Notice 
that matrices C k given by the second relation in (12.21) are positive definite. 
Vice versa, if matrices C k are positive definite, Szego recurrence is uniquely 
recovered from system (11. 2ft . The same is true for the block Szego recurrences 
( 12.11) . where R k > and 9k = 9(R k ) for some function 9. 

Proposition 2.1 There is a one to one correspondence between the subclass 
of system U.ty) , where the matrices Ck > satisfy M.3\) , and block Szego 
recurrences 112. where Rk > and 9k = 9(Rk). This correspondence is 
given by \2.2\) . A2.3\) to map block Szego recurrences into Dirac type systems 
and is given by the equality 

u = i m , R k = (y;c k u k y > 0, u k = u^Rk-i) (2.5) 

to map Dirac type systems into block Szego recurrences. 

Proof. The first part of the proposition is already proved. Moreover, ac- 
cording to (12.21) and (12.31) matrices Rk (k > 0) and Uk (k > 0) are uniquely 



W k (\) 



(i - A-) 



It' 9 



■Uk 



zl p 



L p _ 



Rk 



1 -Pk 
-~P~k 1 
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defined by the relations (12.51) . Clearly we have R k = R k > 0, and it remains 
to prove R k jR k = j- We shall use for this purpose a unitary equivalence of 
U k C k U k to diagonal matrix Dk > 0: 

U* k C k U k = U* k D k U k , U* k U k = U k U* k = I m . (2.6) 

When k > we assume that R r jR r = j is already proved for r < k, and so 
UkjUl = j. Then formula ([21]) implies that UiD k U k iUlD k Uk = j, i.e., 

D k l = JkDkJk, Jk := UkjUl J k = J* k = Jk 1 - (2-7) 

By (12.71) . without loss of generality we can assume that the diagonal matrix 
Dk has the following form: 

D k = diag^o^ dj h , d^I h , d r I lr , d" 1 /^}, (2.8) 

where do = 1 and d\ ^ d s , di ^ d' 1 for I ^ s. Formulae (12.71) and (12.81) imply, 
in their turn, that the matrix Jk has a block diagonal structure 



J k = diag{w , ji, • • • , jr}, js 



u s 
ul 



(2.9) 



where u s are l s x l s unitary matrices, and Uq = u^. According to (12.81) and 
( 12T9D we have 

J k Dl=D k h k . (2.10) 
By (12. 6p we have R k = U k D^U k , and, taking into account (I2.10p . we get 

R k jRk = UiDlJkDlUk = U* k J k U k = U* k U k jU* k U k = j. 



The spectral theory of the discrete and continuous systems is strongly 
related to the construction of the fundamental solutions (see, for instance, 
0-110], ED [MIES], HEQ-H2] and references therein). The j-properties of 
the fundamental solutions play an important role [TQl H21 CEH EES! EH [231 [26J 



For the case of the explicit construction the version of the Backlund- 
Darboux transformation (BDT) introduced in [23 EDI EI] proves very fruitful. 
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Choose n > 0, two n x n parameter matrices A (det A ^ 0) and So = S£, 
and n x m parameter matrix U Q such that 

AS -S A* = M jIl*. (2.11) 

Define sequences {11^} and {S k } {k > 0) by the relations 

n fc+1 = n fe + iA^Ukj, (2.12) 

S k+1 = S k + A~ l S k (A*r l + A-^UKA*)- 1 . (2.13) 
It follows that the matrix identity 

AS k+1 - S k+1 A* = iU k+l3 U* k+1 (fc>0) (2.14) 

is true. Following the lines of the discrete BDT version for the skew self- 
adjoint discrete Dirac type system presented in [21], we get the theorem. 



Theorem 2.2 Suppose det S r ^ (0 < r < N). Then the fundamental 
solution W k+ i of system U.ty) . where 

c k := i m + n^s , A : 1 n fc - ui +1 s^ +1 ii k+1 , (2.15) 

admits representation 

W k+1 (X) = w A (k + 1, A) (J m - lj) fc+1 WA (0, A)" 1 (0 < k < N). (2.16) 
Here W k+ i is normalized by the condition Wq(X) = I m , and 

w A (k, A) := I m - ijU* k 3^\A - Mn)- 1 ^ (2.17) 

The right hand side of ( 12.17ft with fixed k is a so called transfer matrix 
function in Lev Sakhnovich form [35j-[37j. 

We say that system (11.21) . where matrices C k are given by (I2.15p . is 
determined by the parameter matrices A, S and n . 

P r o o f of theorem. Formula (I2.16P easily follows from the basic for this proof 
equality 

WA (k + l,\)(l m -jj) = (L m - l -jC k )w A {k,\)i (2.18) 
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that we shall derive now. Taking into account (12.1 7p one can see that (12.181) 
is equivalent to the equality 

~j(I m - c fc ) =-(i m - l -jC k )ijjr k s^\A - \i n )- l u k 

+ ijUl^S^A - \l n )- l Yi k+1 {l m - (2.19) 
i.e., the Taylor coefficients at infinity of the matrix functions in both sides of 

oo 

(I2.19P coincide. Hence, by the series expansion [A— A/ n ) _1 = — A -1 (A _1 A)' 

r=0 

and formula (I2.12p . formula (I2.19P is equivalent to a family of equalities: 

i m -c k = -n*^- 1 ^ + n* +1 5 fc "in fc+1 (2.20) 

and 

K k A r ~ 2 U k = (r>0), (2.21) 

where 

K k := K +l S k l +1 {A 2 + Q - HIS, 1 A 2 + iC k jUlS^A. (2.22) 

Notice that (I2.20p is immediate from (12.151) . If we prove also K k = 0, then 
(I2.2ip will follow, and so we will get (I2.19P or equivalently (I2.18p . which 
implies (I2.16p . It remains to show that K k = 0. For this purpose we shall 
rewrite fl2T22D using fl2~T2|) and fl2TT5D : 

K k = ni^M 2 + Q - His, 1 a 2 + ijUis^A + iW k s^u kJ w k s^A 

- fll^S^ (n fc + tA~ l U kJ )jUlS^A. (2.23) 

According to (gZEJ) we have ^jll^ 1 = A - S k A*S^ x . Therefore, from 
(I2.23P we derive 

K k = Ul^S^In + S^S^A+A-^UlS^A) -H* k A* S^A+ijU^A. 
In view of (I2.13P we simplify our last formula: 

K k = TT k+1 A*S^A - U* k A*S^A + miS; 1 A. (2.24) 
Finally, by fl2~T2l) and fl2T24j) we have K k = 0. ■ 



Proposition 2.3 Suppose det S r ^ (0 < r < N). Then the matrices C k 
(0 < k < N) given by A2.15\) satisfy conditions 

Proof. The first equality in (11.31) is immediate. To prove the second equal- 
ity notice that by the standard in the S'-node theory [35]- [37] calculations 
(see also, for instance, formula (2.10) in [TTJ it follows from (12. lip and (12.131) 
that 

w A (r,X)*jw A (r,X) = j + i(\ - A)n*(A* - ~\I n y 1 S~ 1 (A - A/ n ) _1 n r . (2.25) 
In particular, we have 

w A {r, X)*jw A (r, A) = j. r> 0. (2.26) 
It is easily checked also that 

( /m+ ^.( /ro _^) = ( 1 + J_),, (2 . 27) 

According to (12.181) formulas (12.261) and (I2.27P yield the equality 

(l m + \c k j)j (l m - % -jC k ) = (l + -I) j. (2.28) 
Therefore the second equality in ( 11. 3ft holds. ■ 

3 Auxiliary propositions 

Recall that the invertibility of matrices S k is essential for our constructions. 
On the other hand the important subcase of Szego recursion corresponds to 
system (11.21) . where C k > 0. A natural condition, when all > and C k > 
is given in our next proposition. 

Proposition 3.1 Let the parameter matrix S be positive definite, i.e., S > 
0. Then we have 

S k >0 (k > 0), C k > (k > 0). (3.1) 
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Proof. The inequalities for Sf~ in (13.11) follow from (I2.13P by induction. To 
derive the relations Cf. > 0, introduce first two block matrices: 



G 



A' 1 aA- l H k 
-ibj 



where 



(3.2) 

a(2 + ac) = l, b{l + ac) = l, (3.3) 

and, moreover, c is sufficiently large so that G > 0. We shall discuss the 
choice of a and b satisfying (13. 3p later on. According to fl 2 . 1 2 j) . (12.131) . (13.21) 
and (13.31) . direct calculations show that 



G + FGF* 



Sk+i 

n 



n 



fc+i 



fc+l c(l + b 2 )I n 



(3.4) 



As G + FGF* > G > we have G' 1 > (G + FGF*)-\ and, therefore, 
the inequality holds also for the m x m right lower blocks of these matrices: 
(G _1 ) 22 > ((G + FGF*)- 1 )^. Finally, we obtain 

(3.5) 



(3.6) 



( G ") 22 j <{((G + FGFT 1 ) 22/ 
Taking into account ( 13.21) . we can rewrite (13.51) in the form 

d m - rr^rL, < c (i + b 2 )i m - ni^s^n^. 

Let us fix c and choose a root a (0 < a < 1/2), of the equation 



2 2 1 
a H — a 

c c 



0. 



which is always possible. Putting also b = a(l — a) -1 , we see that relations 
(13.31) hold. Moreover, the first relation in (13.31) means that a 2 c = 1 — 2a, and 

so 

cb 2 = ca 2 (l - a) 2 = (1 - 2a)(l - a) 2 < 1. (3.7) 
From ([31D and ([HD it follows that 

i m + n^n, - n^s^n^ > o (3.8) 

Recall the definition (12 . 1 5[) to see that inequality (I3.8P implies C k > 0. ■ 
In this section we shall need as well another property of Ct- 
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Proposition 3.2 Let relations U.3\) hold, and assume that C k > 0. Then 
we have C k ± j > 0. 

Proof. It follows from (11.31) that 

(C k + ej)j(C k + ej) = 2e(C k + ^fj)- (3.9) 

If {C k + ej)f = 0, then by igM) we have also (C k + ^fj)f = 0, and so 
(1 — e 2 )jf = 0. Therefore, we have det(Cfc + ej) ^ 0, when |e| < 1. Thus, 
the inequality C k > yields (Ck + ej) > for |e| < 1. ■ 



4 Weyl functions, direct and inverse 

problem: the case of the pseudoexponen- 
tial potentials 

Following definitions of the Weyl functions for Sturm-Liouville, Dirac type 
and canonical systems on the semiaxis (see, for instance, ^Z5\ [37] and refer- 
ences therein), we can define also Weyl functions for system (11.21) . Namely, 
let matrices C k > satisfy (11.31) . Then, a p x p matrix function (p holomor- 
phic in the lower halfplane C_ is said to be a Weyl function for system (11.21) 
on the semiaxis k > 0, if the inequality 



5^b(A)* i p ]qw k Wk(xyCkWk(\) 

k=0 

,2|/|\2| , 



< oo (4.1) 



holds, where q(X) = |A 2 |(|A 2 | + 1) 
Remark 4.1 Similar to the continuous case we have a summation formula: 

J2q(\) k Wk(xrc k Wk(X) = ^±lUxy+^w r+1 (x)*jW r+1 (x)-j). 

(4.2) 

Indeed, according to U.S\) and U.3\) we have 

w k+1 (\)*jw k+1 (\) = w k (\y(i m + =c k j)j(i m - l -jCk)w k (\) 
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q (x)- l w k {xyjw k {x) + \ x) w k (\yc k w k (x), 



i.e., 

|A 2 | + 1 



q(\) h+1 W k+1 (\y 3 W k+1 (\) - q{\) k W k {\YjW k {\) 



i(X - A) 

= q(\) k W k (\yC k W k (\). (4.3) 
Formula h4-3\ ) yields {^.ty . 

To construct the Weyl function, partition first matrix Ilo and matrix-function 
Wa{0, A) into blocks: 



n = [$ w A {o,x) 



a(X) 6(A) 
c(A) d(X) 



(4.4) 



Similar to the considerations in [2JJ it follows from (12.171) that 

b{X)d{\)- 1 = -t<t>*S \A x - XQ- 1 ^, A x = A + iW*S \ (4.5) 
To calculate c?(A) _1 here, we use the fact from the system theory: 

(l p + C(XI n - Ay 1 By 1 = I p -C{XI n -(A-BC)Y 1 B. (4.6) 



Theorem 4.2 Let parameter matrices be fixed, assume Sq > 0, and define 
C k by Ii2.15\) . Then system U.ty) is well-defined on the semiaxis and its 
unique Weyl function, which satisfies \4-l\j , takes the form 

y?(A) = -i<S>*S \A x - XQ- 1 ^, A x = A + iW*^ 1 . (4.7) 

Proof. By Proposition 13.11 system (II. 2p is well-defined. Now, relations (14.51) 
imply if = bd~ l for the matrix function ip given by (14. 7p . According to (12.251) 
we have 

w A (0,xyjw A (0,X)<j (AeC_), 

and it follows, in particular, that d(X)*d(X) > I p + b(X)*b(X). Therefore, we 
get 

p(A)V(A)<ip (AeC_), (4.8) 
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and so ip is holomorphic in C_. Notice that the equality ip = hd 1 is equiva- 
lent to the formula 



<p{\) 



w A (0,X) 







d(X)- 1 . 



(4.9) 



Taking into account (14. 9 p and w^(r + 1, X)*jwA(r + 1, A) < j, we derive from 
representation (12.161) of W r+ i(\) that 



[<p(\)* Ip] W r+1 (A) *jW r+ i (A) 



<p(A) 



|A + z|^ +2 |A 



2r+2| \ \-2r-2 



(d(\y) 



-i 



x [0 I p }w A (r + 1, A)*jw A (r + 1, A) 
By Q£2J and (OOj) the inequality 

r 

b(A)* J p ]^g(A) fe ^ fc (A)*C fc ^(A) 

fc=0 



rf(A)" 1 < 0. 



(4.10) 



<p{\) 



< 



|A 2 | + 1 



i(A- A) 



(4.11) 



is true. From (14. lip inequality (14.11) is immediate, i.e., <p defined by (14.71) is 
a Weyl function. 

Let us show that <p is a unique Weyl function. First notice that by 
Proposition 13.21 we have inequality W*C a W s > W*jW s . Now, use relation 
(14.31) to derive inequality q s W*jW s > q s ~ 1 W*_ 1 jW s -i. From the inequalities 
above we get 

q(\) k W k {\)*C k W k (X)>j. (4.12) 



Therefore, the following equality is immediate for any / e C p : 



J2fV P o]q(\) k w k (\yc k w k (\) 



k=0 







/ = 00. 



(4.13) 



According to (14. ip and (14.131) . the dimension of the subspace L e C m , such 
that for all h G L we have 



J2h*q(\) k W k (\yC k W k (\)h < oo, 



(4.14) 



k=0 
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equals p. Now, suppose that there is a Weyl function ip ^ ip, where is 



given by (14.71) . Then the columns of 



and the columns of 



£(A) 



belong to L. Therefore, dim L > p for those A, where y?(A) 7^ y?(A), an d we 
come to a contradiction. ■ 



Remark 4.3 If Sq > 0, then by Ii2.12\) - I[2.15\) we can substitute parameter 
matrices A, Sq and Uq by the parameter matrices S 2 AS 2 , I n and S 2 Ho, 
which determine the same system. For Sq = I p formula \4-ty takes the form 



XL 



A x = A + i^^*, 



(4.15) 



and we have also A x - (A x )* = i($$* + tftf*), det(A x - i^m*) ^ 0. 



Example 4.4 Consider the simplest example: p = 1, n = 1, A = a & 
(a ^ 0), S = 1. From KETT\) and it follows that |$| = |^| and 



IT 



a + i 



% U k W k = 2|$| 



a 2 + 1 \ * 



(4.16) 



iVoit;, m meu> 0/ So = 1, \2.1$) and the second relation in ( [^.i6] j one can 
check that 



c 



a 2 + l' 



(4.17) 



Finally, using /12.15\) . {4-10$ and \4-l r /\) we oe£ t/ie entries (Cfc)y 0/ CV 

(C fe )n = (C,) 22 = 1 + C|* | 2 (fcC + ^"'((^ + !)C + (4-18) 



'a + i 



fc+i 



(Cfcjai = (c,) 12 = m ((*C + l)" 1 - ((* + 1)C + 1)- 1 ( 

(4.19) 

T/ie Weyi function of system ( f i ■-§)) ■ where the matrices Ck are given by H4-18$ 
and \4-iyj) > i s easily calculated using fl^.i5] ): 



(/3(A) = z<M>(A - a - 



2\-l 



(4.20) 



Notice that our matrices Ck belong to the class of the so called pseudoex- 
ponential potentials. An important subclass of the strictly pseudoexponential 
potentials, that is, a subclass with an additional requirement o~(A) C C_ (a - 
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spectrum), have been treated for p = 1 in [5JE]- In particular, for the strictly 
pseudoexponential subcase the inequality |y?(A)| < 1 for A G C_ is true. On 
the other hand, in the simple example above we have cr(A) = a G R and 
\(p\ — 1 for A = a. 

According to ( 14. 7ft and fl4.8D Weyl function ip is a rational, strictly proper 
and contractive in C_ matrix function. By the proof of Theorem 9.4 [22] such 
matrix functions admit representation (realisation): 

cp(X) = —i$*(9 - A/ n ) -1 ^, (4.21) 

where 

0-9* = i($$* + (4.22) 

Direct calculation shows also that formulas (14.21 j) and (14.221) yield I p — ip*(p > 
for A G C_. So, realization ( 14.2 ip . (14.221) is equivalent to function being 
rational, strictly proper and contractive in C_. 

Theorem 4.5 Matrix function <p is the Weyl function of some system U.2\) 
determined by the parameter matrices with Sq > if and only if it admits 



representation H4-21\) , d4-2ty such that det(8 - i^^*) ^ 0. Then y? is the 
Weyl function of some system U.2\) . where Ck > 0. To recover such system 
put 

S = I n , A = 9-im*, $ = $, * = n = [$ (4.23) 

and define matrices Ck by formula Ii2.15\) . where matrices Uk and Sk (k > 0) 
are given by formulas A2.12]) and A2.13\) . 

Proof. The necessity of theorem's conditions follows from Remark 14.31 
Now, suppose these conditions are fulfilled. Then, from (I4.22j) and (14.231) it 
follows that the identity (12.111) holds for the parameter matrices. Therefore, 
system ( ll.2p is defined. So, by Theorem 14.21 ip is the Weyl function of this 
system. ■ 

Remark 4.6 The Weyl functions in the upper halfplane can be treated in a 
quite similar way. That is, we define Weyl functions in C + by the inequality 



J2[i p ^\y}q(\) k w k (\yc k w k (\) 



k=0 



MA) 



< oo. (4.24) 
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Then the Weyl function of system II 1.2$ . where matrices C k are given by 
112.15}) and S > 0, takes the form 

(p{\) = c(A)a(A)- 1 = i$*S' - 1 (A x - XIn)- 1 ®, A x = A- i<5><5>*S \ (4.25) 

A definition of a Weyl function in C_ can be also given in a more general 
form. 

Definition 4.7 Let matrices C k > satisfy U.3\) . Then, a p x p matrix 
function holomorphic in C_ is said to be a Weyl function for system M.2\) 
on the semiaxis k > 0, if the following inequality holds: 



k=0 



-icp(X) 



< oo. 



(4.26) 



Here K* = K~ l and q(X) = |A 2 |(|A 2 | + 1) 



When K in the inequality (14.261) equals I n , this inequality coincides with 
the inequality (14.11) . In general, the choice of the matrix K is related to 
the choice of the domain of the operator corresponding to the Dirac system, 
and usually K is chosen so that the Weyl functions are Herglotz functions. 
Further we assume that 



K 



1 

7! 



ip ip 



(4.27) 



Simple transformations show that the Weyl function ifj defined via (14. ip and 
the Weyl function (fx defined via (14.261) and (14.271) are connected by the 
relation 

<p K = -t(I p -<p I )(I p + <p I )- 1 . (4.28) 
From ffl~2gj) it follows that 

MA)-MA)* = -2t(I p + Vl (\)T\lp-¥>iW*ViW)(Ip+ViW)~\ (4-29) 

A G C_. Thus, according to H4.8[) and (14.291) fx is a Herglotz function with 
a non-positive imaginary part in C_. 
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5 Weyl functions, direct and inverse problem 
on the interval: general case 

In this section we shall consider the self-adjoint matrix discrete Dirac type 
system (11.21) on the interval < k < N. We assume that (11.31) holds and C k > 
0. It was shown in Proposition [2J] that C k > yields C k = (U k *)~ 1 RlU k 1 , 
where R k = R* k and R k jR k = j- Hence, we get 

C k = (Uty\Rl + R k jR k )U k l — j = (U* k )- l R k (I m +j)R k U k l -j 

= 2p(k)*p(k) - j, < k < N, 

where 

m = [I p 0}R k U k \ P(h)jP(k)* = I p . (5.1) 
Further we shall use these relations: 

C k = 2/3(k)*P(k) - j, P(k)j/3(ky =I P , < k < N. (5.2) 



Remark 5.1 Relations 115. fy) are equivalent to the relations C k > and 
U.3\) . Indeed, we have just derived $5.0j) from C k > and C k jC k = j , and 
vice versa: direct calculation shows that A5. yields To derive also 

from ( 15. S\) the inequality C k > 0, choose a matrix (3{k) such that 



P(k)jp(ky = o, p(k)jp(ky 



-I,. 



(5.3) 



Notice, that in view of the second relations in A5.2\) the maximal subspace, 
which is j- orthogonal to the rows of (3 k , proves to be p- dimensional and j- 
negative, i.e., (3(k) always exists. According to 115. fy) and ( 15. ty) we have 



p(k) 



Pik) 
pik) 



pjk) 
p(k) 



Pik) 
P(k) 



(5.4) 



Finally, by the first relations in K5.2§ and by \5.1$ we obtain 



c k = p(kyp{k) + p(kyp{k) 



' P(k) ' 




" P(k) ' 


_ P(k) _ 




P(k) 



> 0. 



(5.5) 
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From the second relation in (15. 2p for k > it follows also that 

det (j3(k)jfi(k + 1)*) ^ 0, < k < N - 1. (5.6) 

Indeed, if (15.61) does not hold, we have (3(k)jf3(k + 1)*/ = for some f^O. 
Then, in view of the second relations in (15.21) for k > 0, we see that the linear 
span of the rows of (3 k and of f*/3k+i forms ap + 1-dimensional j-positive 
subspace of C m , which is impossible. 

Notice also that according to (12. ip and (12.31) we have 



-ijR^U^, U* k = -iR^U^j (k > 0). 



(5.7) 



Compare (15.11) and (I5.2p to see that a matrix (3{k — 1), which satisfies the 
equalities 



C, 



k-l 



2p(k - iyp(k - 1) - j, p(k - i)jp(k - iy 



coincides with the upper block row of —iRk-iU^li up to a p x p unitary 
factor u(k — 1). Moreover, (3(k — 1) given by (I5.3P defines the lower block 
row of —iRk-\U k \ up to some p x p unitary factor u(k — 1). Taking into 
account the second equalities in (12. 5p and (15.71) . we get 



R k = \^D{k-l) 
where 



l(k - 1) 
(3(k-l) 



D(k-1) 



jc k3 p(k-iy p{k-iy D(k-i) 



u{k - 1) 
u(k-l) 



(5.8) 
(5.9) 



Remark 5.2 We see that the two matrices C k and C k -\ determine R k by 
Ii5.8\) up to a unitary block diagonal matrix D(k — 1). 

Similar to the continuous case, Weyl functions of the discrete system on 
the interval are defined via Mobius (linear-fractional) transformation 

cp(X) =z(W 21 (A) J R(A)+W 22 (A)g(A))(W 11 (A) J R(A)+W 12 (A)g(A))^ 1 , (5.10) 
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where R and Q are pxp analytical functions in the neighbourhood of A = — i, 
and 

W(A) = {Wij(X)}l j=1 = KW N+1 (X)\ (5.11) 

Here, coefficients WV/ of the Mobius transformation are the pxp blocks of 
W, the matrix K is given by (I4.27P and 



K* = K-\ KjK* = J, J 



I p 
I P 



(5.12) 



It would be convenient to put (3{k) := j3{k)K* and rewrite (15. 2p as 



C k = 2K*P(k)*P(k)K - j, P(k)jp(k)* 



< k < N. 



(5.13) 



We shall need the following analog (for the self-adjoint case) of Theorem 3.4 



Theorem 5.3 Suppose W (Wq(A) = I m ) is the fundamental solution of 
system U.ty) , which satisfies conditions A5.13\) . Suppose also that a p x p 
matrix function ip is given by formulas li5.1U\) and Ii5.11\) . where 

det(VM-i)i2H) + Wi 3 (-i)Q(-i)) 0- (5.14) 
Then system U.ty) satisfies U.HJ\) , Ck > (k >0), and the inequalities 

det (j3(k)Jf3(k + 1)*) ^ 0, < k < N - 1 (5.15) 

hold. Moreover, system U.2fy is uniquely recovered from the first N + 1 Taylor 



coefficients {ak}k=o of icp I i 



z + 1 
z- 1 



at z = by the following procedure. 



First, introduce (N + l)p x p matrices $i, $2 : 



a 









$1 = 




, $2 = 









a + ai + . . . + on 



(5.16) 
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Then, introduce an (N + l)p x 2p matrix U and an (N + l)p x (iV+ l)p block 
lower triangular matrix A by the blocks : II = [$! $ 2 ], 







for r > 



A:=A(JV) = {a 3 _ k } 



N 

kJ=o 



-I p for r = . (5.17) 
i I p for r < 



Next, we recover (N + l)p x (TV + l)p matrix S as a unique solution of the 
matrix identity 

AS - SA* = mm*. (5.18) 

This solution is invertible and positive, i.e., S > 0. Finally, matrices (3{k)*f3{k) 
are easily recovered from the formula 



WS~ l Ii = B*B, B := B(N) 



0(0) 
0(1) 

(3(N) 



(5.19) 



/Vow, matrices Ck and system are defined via the first equality in Ii5.13\) . 

Proof. Step 1. According to Remark 15. II relations Ck > and (11. 3p follow 
from (I5.13p . Relations (I5.15P follow from (15.61) . Now, let 



K(r) 



Mr) 
K x {t) 

KJr) 



(5.20) 



where Ki(r) are p x (r + l)p matrices of the form 

K l (r) = i/3{l)J\j3(0)*.../3(l-iy P(l)*/2 0...0]. (5.21) 
From d5TT9|) - flCT]) it follows that 

K(r)- K(r)* = iB(r)JB(r)* . (5.22) 
By induction we shall show in the next step that K is similar to A: 

K(r) = y_(r)v4(r)V_(r) _1 (0 < r < TV), (5.23) 
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where V_(r) ±:L are block lower triangular matrices. Taking into account 
(I5.23P and multiplying both sides of (I5.22p by V_(r) _1 from the left and 
by (VL(r)*) from the right, we get 

A(r)S(r) - S(r)A(r)* = iU(r)JU(r)*, (5.24) 

S(r) := V_(r) _1 (VL(r)*)~\ U(r) := 7_(r) -1 B(r). (5.25) 

Moreover, Step 3 will show that matrix V-(N) can be chosen so that the 
equality 

n = [$! $ 2 ] = V^(N)~ 1 B(N) (5.26) 

holds, i.e., II = II(iV). (Here <£>i and $ 2 are given by (15.161) .) 

Identities ( I5.24p have unique solutions S(r) as the spectra of A(r) and 
A(r)* do not intersect. (The statement follows from the rewriting of (15.241) 
in the form 

S(r)(A(r)* - XI)- 1 - (A(r) - Xiy^r) 

= %{A{r) - AJ) _1 n(r) m(r)*(A(r)* - XI)- 1 , 

and from the following integration of both sides of the obtained identity along 
a contour, such that the spectra of A is inside and the spectra of A* outside 
it. In particular, by (15.181) and (I5.24p one can see that S := S(N). Hence, 
we derive from (15.251) and (15.261) that S > and the first equality in (15.191) 
holds. It remains only to prove (15.231) and (I5.26p . 

Step 2. Now, we shall consider block lower triangular matrices V~(k) 
(0 < k < N): 



VC(0)=v_(0) = ) 9 1 (0), V-{k) 



V-{k-l) 
X(k) v.(k) 



(k > 0), (5.27) 



where v^(k) are p x p matrices, where (3\{k) and ^(k) are p x p blocks of 
(3{k) = \fti(k) foik)}, and where X{k) = [X (k) X(k)] are px kp matrices. 
Here Xo(k) are arbitrary pxp blocks, and the matrices X(k), V-(k) are given 
by the formulas 



X(k) = i( (3(k)J[(3(0)* . . . (3(k - l)*]V^(k - 1) 



'(jfc-l)p 




- v_(k)[I p . . . I p 
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x \ A(k - 2) + -/(fc-ijpj , v-(k) = (3(k)J(3(k - l)*v-(h - 1). (5.28) 

According to (15.171) we have A(0) = (i/2)I p . From the second relation in 
( E3SJ) and definitions (ICTj) and <KT^ it is immediate that K(0) = (i/2)I p , 
and so (15.231) is valid for r = 0. Assume that (15.231) is true for r = k — 1, and 
let us show that (15.231) is true for r = k too. It is easy to see that 



V-(k)- 1 



V.(k-l)- 1 
_ —V-(k)~ 1 X(k)V-(k — l)" 1 v-(k)~ 1 



(5.29) 



Then, in view of definitions (I5.17P and (15.271) . our assumption implies 



V-(k)A(k)V-(k)~ 1 



K(k 



Y(k) 





i-I 

2 P 



(5.30) 



where Y(k) = (X(k)A(k - 1) + iv_(k)[I p . . . I p }) iu_(jfe) 

VL(fc-l)- 1 
X _ -u_(A;) -1 X(A;)V r _(A; - l) -1 _ ' 

Rewrite the product on the right-hand side of the last formula as 

Y{k) = (x{k)(A{k - 1) - % -I kp ) + iv.{k)[I p ...I p ^V-[k - 1)- 
From (15TT7j) and fl5UT|) it follows that 



(5.31) 



Y(k) = ^x(k)(A(k-2)+~i {k „ 1)p )+iv4k)[i p ...i p }) ^(^Jy-^-i)- 1 . 

(5.32) 

Notice that the sequence [I p . . . I p ] of identity matrices in (15.321) is one block 
smaller than in (I5.3ip . By (15.281) and f!5.32j) we have 



Y(k) =i(3{k)J [p(0f . . . P(k-l)*)V-(k-l) 



l (k-i)p 




Finally, formulas (15.271) and (15.331) imply 

Y(k) = ip(k)J\p(0)* ■ ../3(k - 1)*] (k > 0). 



f3{k-l)*v_(k-l) 
(5.33) 



(5.34) 
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According to the second relation in f 1 5 . 1 3 j) and formulas (15.211) and (I5.34p we 

get 



Y(k) -I P 



K k {k). 



(5.35) 



Now, using f)5.20p and (15.351) . one can see that the right-hand side of (15.301) 
equals K{k). Thus, ( 15.231) is true for r = k and, therefore, it is true for all 
< r < n. 

Step 3. To derive (I5.26P we shall first prove that the matrices V_(r) given 
by (I5.27P and (15.281) can be chosen so that 



V^r^B^r) 



Bi{r) := B{r) 







Pi(0) 
I Mr) J 



(5.36) 



In other words, the blocks Xo(r), arbitrary till now, can be chosen so. Indeed, 
by the definition in (I5.19P and the first equality in (I5.27P formula (15.361) is 
true for r = 0. Assume that (I5.36P is true for r — k — 1. Then, from (15.291) 
it follows that (15.361) is true for r = k, if only 



v_(k)- l X{k) 



+ v4k)- 1 (3 1 (k) = I p . 



(5.37) 



It implies that we get equality (I5.36P for r = k by letting 



X (k) = Px(k)-v-(k)-X(k) 



(5.38) 



Hence, by a proper choice of the matrices A (r) we obtain (I5.36P for all 
r < N. 

It remains to prove that 



V^(N)- 1 B 2 {N) = $ 2 , B 2 {N) :-- 



&(0) 



(5.39) 



For that purpose we shall consider the matrix function Wjv+i(A), which 
is used in ( 15. lip to define the coefficients of the Mobius transformation 
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(I5.10p . Namely, we shall prove the transfer matrix function representation of 

w N+1 (\y. 

^Pj K*w A (N,-^)K, (5.40) 

where 

w A (r, A) = I 2p - im{ryS{r)- l (A{r) - AJ^)"^). (5.41) 
Identity (15.241) implies a similar to (I2.25P equality 

w A (r,/jyjw A (r,\) = J 

+ i(£-A)n(r)*(A(r)*-/ZJ (r+ ^ (5.42) 

Moreover, according to factorization theorem 4 from [35] (see also [37], p. 
188) we have 

w A (r,X) = (l 2p -im(ryS{r)- l P*{PA(r)P* - \I p y l {PS^P*)' 1 

x PS{r)- l Tl{r)\w A {r - 1,A), P = [0 ... ZJ. (5.43) 
Taking into account (I5.17p . (I5.25p . and (15.271) we obtain 

(PA(r)P* - XIp)' 1 = (^-A)- 1 ^, P£(r)- x P* = v_(r)*v_(r), (5.44) 

P5(r) _1 n(r) = u_(r)*PP(r) = v_(r)*/3(r). (5.45) 
Substitute f[5T34"P and (ET35j) into (ET35I) to get 

u*(r, ^) = (lap - ^JW(r))ti; A (r - 1, (5.46) 
From the definitions (I5.17p . (15.251) . and (I5.4ip we also easily derive 

w A (0, ~) = Jap - Jl- JP(0)*P(0) = J 2p - JL- J^(0)*/3(0). (5.47) 
2 2 — A « — A 

On the other hand system (11.21) with additional conditions (15.131) can be 
rewritten as 

W(r + 1, A) = ^(h P - ^jK*(3(ryp(r) K) W(r, A). (5.48) 
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In view of the normalization W(0) = l2 P , formulas (I5.46I) - (I5.48I) imply (l5.4Up . 
From (ET3Q|) and <^M> it follows that 



W(N + l,X)jW(N+l,X)* 



A + i 
A 



N+l 



X-i 
X 



N+l 



J- 



(5.49) 



Let us include functions <p into consideration. Introduce 

2N+2 



A(X) := 



A 



A + i 



MA)* J p ]iirW(JV+l, A)*jW(iV+l, A)iT 



-^(A) 



(5.50) 



According to floTTUl) . f[5TTTj) . and f[5T4l?j) we have 



A-i 



A 



2N+2 



{(W n (X)R(X) + W 12 (X)Q(X)Y) 



*\ -1 



x (i2(A)*i2(A)-g(A)*Q(A))(Wu(A)i2(A)+Wi 2 (A)Q(A)) \ (5.51) 
By ( 15.14[) and ( 15.511) ^4 is bounded in the neighborhood of A = — i: 

\\A{X)\\ = 0(1) for X^-i. (5.52) 

Now, substitute (ET30j) and into fl530l to obtain 

.4(A) = MA)' ^(j+^A-A)^^)*^)^^/^)- 1 ^)- 1 



x (A(iv) + ^ (iV+1)p )- 1 n(iv) 



-iy?(A) 



(5.53) 



Notice that S(N) > 0. Hence, formulas (15321) and (153B1 imply that 



(A(iv) + ^/ {JV+1)p )- 1 n(iv) 



-z</?(A) 



0(1) for X^-i. (5.54) 



Recall that II(iV) = V_(A^)~ 1 J B(A^) and that A(JV) is denoted by A. Now, 
represent n(iV) in the block form 



I1(N) = [$ 1 (N) $ 2 (iV)], ® k (N) = V4N)- 1 B k (N) (A; = 1,2). (5.55) 
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According to (15.161) and (15.361) we have $i(iV) = $1. Hence, multiplying the 
matrix function on the left-hand side of (I5.54j) by i ( (A+±I {N+1)p ) ~ $ x J $^ 



we derive 



for A 



(5.56) 



The matrix A+^I^ + i) p is easily inverted explicitly (see, for instance, formula 
(1.10) in [32]). As a result one obtains 



$*(A + ^J (7V+1)P )" 
Moreover, we get 



% + A 



[q 



pjj 



A-z 

A + ^ V 
(5.57) 



A 2 

$t(A + -v +1)p )- 1 $ 1 = - TI 



if + 1 - J. 



-i 



(5.58) 



Let A 



z + 1 
^ - 1 



i.e., z 



A + i 
A-z 



. Then, we derive from (15.581) that 



(^(^+^V + i)p) _1 ^i) _1 =(-^ +1 + 0(^ +2 ))/ p (*-0). (5.59) 



Taking into account (I5.57P and (15.591) . we rewrite (I5.56P as 



(z(^))+i(l-z)[/ p *I, z 2 I p ...]$ 2 (iV) 



0(2 



(5.60) 



for z -»• 0. From (IBTTHD and (IOTP it follows that $ 2 (JV) = $ 2 , i.e., (jESSj) is 
true. As $i(iV) = $ a and $ 2 (iV) = $ 2 , so II(iV) = II and formula (l5T26|) is 
finally proved. ■ 

Definition 5.4 Let matrices Cj, satisfy A5.13\) . Then, apxp matrix function 
ip holomorphic in C_ is said to be a Weyl function for system on the 

interval < k < N, if (p admits representation $5.101) . where a pair R, Q is 
meromorphic in C_ ; well-defined at A = —i, and nonsingular with j -property , 
i.e., 



R{\YR{\) + Q(A)*Q(A) > 0, R{\YR{\) < Q{X)*Q{\). 
The set of Weyl functions is denoted by J\f(N). 



(5.61) 
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Using notation (15.111) . we deduce from (14. 2 p the inequality 

q(\) N+1 W(\)jW(\)* < J, A G C_. 
According to [26] we can change the order of factors in (15.621) : 

q(\) N+1 w(xyjw(\) < j, ag c_. 

Moreover, after excluding A = — i the inequality is strict 

q{\) N+1 W{\y JW(X) < j, Ag C_\ - i. 
In view of (Oil . lETTTj) . (IQjl . and ([53]) at A = -i we get 



W(-i) = KW N+1 (z)* = {-2) N+1 Kj[(ftP, 



k3, 



(5.62) 



(5.63) 



(5.64) 



(5.65) 



fc=0 



From the second relations in ( 15.31) and from ( 15.611) we, analogously to the 
proof of (15.61) . derive: 



det[/ p IJj^fl, det P k j0* k+1 ± 0, det(3 N j 
By (15.651) and (15.661) the next proposition is valid. 



R(-i) 

L QH) J 



+ 0. (5.66) 



Proposition 5.5 Let the pair R, Q satisfy ( 15. 61\) Then inequality \5.1J$ is 
fulfilled. 

By Proposition 15.51 and the proof of Theorem 15.31 we get a corollary. 

Corollary 5.6 Weyl functions of system Al.fy) , which satisfies conditions 
h5.13\) , are Herglotz functions and admit the Taylor representation 

K*(l3T)) = - i (^+(^i-i'o)z+. . .+(^n-^n-i)z n )+0(z n+1 ), (5.67) 



where z — > and the p x p matrices ipk ore the blocks of 

<I>., = ' = V-(N)~ 1 B(N). 



N 



(5.68) 
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Proof. From (15301) . (l53Tj) and (I5^3j) it follows that 



I, 



p 



<o, 



— zip 



i.e., 9<^(A) < for A G C_, and so ip is a Herglotz function. 

By Proposition 15.51 the Weyl functions satisfy conditions of Theorem 15.31 
Then, by the second relation in (15.161) we have representation (I5.67P of ip via 
the blocks of $2- By the proof of Theorem 15.31 we get also $2 = ^2(^)5 i-e., 
( I5.68P holds. Here V-(N) and B(N) are recovered from the matrices (3{k) 
and do not depend on (p. ■ 

Remark 5.7 As Weyl functions ip satisfy conditions of Theorem \5.3[ so the 
procedure given in Theorem 15.31 provides a recovery of system U.S\) from a 
Weyl function (i.e., provides a solution of the inverse problem). 

The following proposition is also true 

Proposition 5.8 The set Af(N) (N > M) is imbedded in N(M), i.e., 
N{N) C N(M). 

Proof. By (I4.2p we have 



q(\f +1 W N+1 (\yjW N+1 (\) < q(\) M+1 W M+1 (\yjW M+ i(\), A G c 



(5.69) 



Insert the length N of the interval into the notation W: 



W(N, A) = W(A) = KW N+1 (\)*. 



(5.70) 



From (15T69D and (IBTTOl) it follows that 



q(\f- M (W(M, X)- l W(N, A)) *jW(M, X)- l W(N, A) < 3 



(5.71) 



Moreover, in view of (11.21) and (15.701) we have 



W(M,\y 1 W(N,\)= f[ (i m +{c k j), 



N 



(5.72) 



k=M+l 
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and the expression on the left-hand side of (15.721) is analytic at A = —%. 
Suppose now that ip G Af(N) is a Weyl function generated by some pair R, 
Q, which satisfies (I5.6ip . Then, according to (I5.6ip . (I5.7ip and (15.721) the 
pair 



R(X) 

QW 



R(X) 
[ Q(A) 



(5.73) 



= W(M, A) _1 W(iV,A) 

satisfies conditions of Definition 15.41 too. Moreover, it is easy to see that 
i(Wai(M, X)R{X) + W 22 (M, A)Q(A)) (W n (M, X)R(X) + W 12 (M, A)Q(A)) _1 , 

= i(W 2 x(N, X)R(X) + W 22 (iV, A)Q(A)) (Wn(N, X)R(X) + W 12 (iV, X)Q(X))~ 1 

(5.74) 

= V(A), 

which completes the proof. ■ 

Theorem 15.31 and Proposition 15.81 imply a Borg-Marchenko type result: 

Theorem 5.9 Let <p and (p be Weyl functions of the two discrete Dirac type 
systems M.2^l , which satisfy conditions Ii5.13\) . Denote by Ck (0 < k < N) the 
potentials Ck of the first system and by (0 < k < N ) the potentials of the 
second system. Denote Taylor coefficients of i<f^i(^^\j and iip(i(^^j} 
at z = by {5^} and [ak], respectively, and assume that = ctk for all 
k <N < min{iV, iV}. Then we have C k = C k for k < N. 



Proof. According to Proposition [578], (p and tp are Weyl functions of the first 
and second systems, respectively, on the interval < k < N. By Theorem 
15.31 these systems on the interval < k < N are uniquely recovered by the 
first N + 1 Taylor coefficients of the Weyl functions. ■ 

An interesting Borg-Marchenko type result for supersymmetric Dirac differ- 
ence operators have been obtained earlier in [IT] . 



28 



6 Toeplitz matrices and Dirac system on the 
semiaxis 



By |28j, p. 116 it is easy to recover a block Toeplitz matrix S which satisfies 
f !5.18p . where the blocks $1 and $2 of II are given by (15.161) . Namely, we 
have 

S = {sj_fe}^ =0 , s_ k = a k = s* k (k>0), s = s* = a + a* Q . (6.1) 

Moreover, this 5* is a unique solution of (I5.18p . A description of all exten- 
sions of S preserving the number of negative eigenvalues, which uses transfer 
matrix function wa, is given in [28] (see also Theorem 4.1 in [32]) in terms 
of the linear fractional transformation 

<p{X) = (R{X)w n (X) + Q(A)w 21 (A)) 1 (r(X)w 12 (X) + Q(A)™ 22 (A)) , (6.2) 

where {wkj{X)} 2 k ^ =l = wa(N,X), and the meromorphic pairs R, Q have J- 
property, i.e., 

R(X)R(X)* + Q(X)Q(X)* > 0, R(X)Q(X)* + Q(X)R(X)* > 0, AGC+. 

(6.3) 

In particular, for the case S > 0, which is treated here, the matrix func- 
^( i(z + 1) \ 

tions tpy — —. -j-J are always analytic at z = and admit the Taylor 

representation 

^( i(z + l)\ ^ ^ ^9 / 
V?( — 7T7 + ) = s + s^z + s- 2 z + ■ ■ ■ (6.4) 



2(z-l 

Our next statement follows from Theorem 4.1 |32|. 



Theorem 6.1 Assume that S = {sj-k}kj=o > ®> an d fi x a o such that 
oto + «q = so. Using 115.16}) , \5.J^1% and 116. 1\) introduce II = [$1 $ 2 ] and 

{wkj{X)}lj =1 = wa(N,X). Now, let matrix functions (p be given by $6. 

where the pairs R, Q satisfy ^6.3\) and are well defined at X = ~. Then 

^ ^ ( i(z + 1) 
the Taylor coefficients at z = of the matrix functions (py — ^ 

satisfy relations 

= s_ fc (0 < k < N), s = a . (6.5) 
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Moreover, putting s_ 



slit for k > N, we have {sj-&}fci=o — /• 



0'/' 



all M > N . In other words, the Taylor coefficients of the matrix functions 
^( i(z + 1) \ . „ . ,, , 

V~ 2~( 1)/ y enera ^ e nonne 9 a ^ we extensions of b. All the nonnegative 

extensions of S are generated in this way. 

Taking into account that wa{N, X)Jwa(N, A)* = J, we derive equality tp = <p 
for the matrix function 



£(A) = -(w l2 (X)*R(X) + w 22 (X)*Q(X)) (wn(\)*R(\) +w 21 (X)*Q(X) 
where 



(6.6) 



R(X)*R(X) +Q(A)*Q(A) > 0, R(X)Q(X) + Q(X)R(X) = 0, A e C_ 



Notice also that relations (16.31) and (16.71) yield 

J R(A)*Q(A)+g(A)* J R(A) <0, 



(6.7) 



(6.8) 



and vice versa relations (16 ,7p and 



yield the second relation in (16. 3p . 



Hence, Theorem 16.11 can be reformulated in terms of the linear fractional 
transformations (16. 6p . where R, Q have J-property (16. 8p . Finally, use (15. lip . 
f l5TT2|) and fl5\i0|) to rewrite floUOl) in the form 



MX) = -(w 12 {-X/2)*R{-X/2)+w 22 {-X/2)*Q{-X/2) 



x ( Wll {-X/2)*R{-X/2) + w 21 {-X/2)*Q{-X/2) 



where we put 



' R{- 


-A/2) " 


= K 


- R(X) - 




-A/2) _ 


. Q(A) . 



(6.9) 
(6.10) 



Here, formula (16.101) is a one to one mapping of the pairs satisfying (15.611) 
into pairs satisfying the first relation in (16. 7p and relation (16. 8p . By (16.61) and 

i(z + l)\ ~( i(z + l)\ . f.(z + l)> 



we have (p(^ — 



Theorem 



2{z-l) 
can be rewritten. 



n - 



2(z-l) 



Therefore 
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Theorem 6.2 Assume that S = {sj-k\kj=o > 0> a o su °h that a + «q = 
So, and introduce W via A5.ll]) and ( |5.^| ). Let matrix functions ip be given 
by A5.10\) . where the pairs R, Q satisfy A5.61\) and are well defined at A = —i. 
Then itp(—i) = ao, and the following Taylor coefficients a^ at z = of the 

matrix functions i(f(i-, — ~~Tj) sa ^ s f v re ^ a ^ ons 

a k = s_ k (0<k<N). (6.11) 

Moreover, putting S-k = s t = a k for k > N, we have {sj-k}kj=o — f or 
all M > N . In other words, the Taylor coefficients of the matrix functions 

^(*(~~~~T)) y enera ^ e nonne 9 a ^ ve extensions of S. All the nonnegative ex- 
tensions of S are generated in this way. 



Remark 6.3 By Definition \5.4\ and Theorem \ 6.S\ the Weyl functions from 



the Weyl diskM{N) generate all the nonnegative extentions of S. It provides, 
in particular, another proof of Proposition 15. 8\ 

(6.12) 

Consider now system (11.21) . which satisfies ( 15. 13ft on the semiaxis k > 0. 

Theorem 6.4 Let system U.fy) be given on the semiaxis k > and let ma- 
trices Ck satisfy Ii5.13\) . Then, there is a unique function which belongs 
to all the Weyl discs J\f(N) : 

oo 

f)fS(N)=< Poo . (6.13) 

N=0 

Proof. According to Corollary 15.61 matrices {Ck}^ =0 (N < oo) or equiva- 
lently matrices {/5(fc)}^ uniquely define blocks {s_fc}^ , where = a k = 
— i>k-i for > and So = «o + «o ( a o — V'o)- Moreover, by Proposi- 
tion 15.81 these S-k do not depend on iV > k, and so system ( 11.21) on the 
semiaxis determines an infinite sequence {s~k}k^=o- By Theorem 15.31 we have 
{ s j-k}kj=o > f° r a ^ N >0. Apply now Theorem 16.21 to see that 



^tjtit) =ad + E s -*^ = < ^( < feriy)' (6 - 14) 
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i.e., this <f belongs to O^^Ni^N). Moreover, as the sequence {s_fc}^ is 
unique, so by Theorem 16.21 the function ip e HivLo -^(-^O * s un iQ ue - ' 



Recall that a Weyl function on the semiaxis is defined by Definition 14. 7\ 
where K is given by formula (14.271) . Theorem 16.41 yields our next result. 

Theorem 6.5 Let system U.2fy be given on the semiaxis k > and let ma- 
trices Ck satisfy A5.13\) . Then, the matrix function y?oo given by A6.13\) is the 
unique Weyl function of system U.ty) on the semiaxis. 



Proof. By (15~TUj) and (ISTTSj) . we have 



-i^oo(A) 



JW(r + l,\) 



R(X) 



(6.15) 



for all r > and for some depending on r pairs R, Q, which satisfy (I5.6ip . 
In view of ( IBTTTjl . (15321) . fl5^9|) . and fl6U5|) we obtain 



[ifooixy i P ] (q(\y +1 Kw r+1 (\y 3 w r+1 (\)K* - j 



-i^oo(A) 



i((Poo{>) - <^oo(A)' 



+ 



| A 2 + 1 |2 X^+l 



|A|W + 1 
Now, formulas (I5.6ip and (16.161) imply 

[ipoo(A)* I P ] (q(X) r+1 KW r+1 (\yjW r+1 (X)K* - J 

< «(<Ax>(A) - ^oo(A)*V 

It follows from (TO]) and flBTTTD that 

r 

X>Voo(A)* / p ]g(A) fc iriy fc (A)*C fc iy fe (A)ir* 



[i?(A)* Q(X)*]j 



R(X) 

Q(A) 
(6.16) 



-i^oo(A) 



(6.17) 



fc=0 



-^oo(A) 



|A| 2 + 1 

< (^oo(A) -^oo(A)*). 



(6.18) 
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Finally, by (16.181) the inequality (I4.26P is immediate, and given by (16.131) 
is a Weyl function. 

To prove the uniqueness of the Weyl function notice that by Proposition 
and by relation (14.31) we have 



q k W* k C k W k > q k Wt 3 W k > q^'W^jW^ >...> W*jW =j (A G C_). 

(6.19) 

Hence, in view of (16.191) we obtain 

J2[I P I p }q^) k KW k {\yC k W k {\)K* \ Ip 1 > 2(r + 1)1 



k=0 

and it follows that 



J2[i P i P ]qW k Kw k (\yc k w k (\)K* 



k=0 



oo. 



(6.20) 



Taking into account Definition 14.71 and inequality (I6.20p . we can show the 
uniqueness of the Weyl function similar to the proof of the uniqueness in 
Theorem 14.21 ■ 



Now, we formulate a solution of the inverse problem. 

Theorem 6.6 The set of the Weyl functions of systems M.Sfy , given 
on the semiaxis k > and such that matrices C k satisfy Ii5.13\) , coincides 
with the set of functions ip such that 



(£+1) 
(z-1) 



aco 



+E 

k=l 



S-kZ 



(6.21) 



are Caratheodory matrix functions in the unit disk and {sj- k } k j =0 > for 
all < iV < oo (s := a + «g). These systems M.ty) are uniquely recovered 
from their Weyl functions via the procedure given in Theorem 15.31 

Proof. According to Theorem 16.51 the Weyl function on the semiaxis is also 
a Weyl function on the intervals. Hence, the procedure to construct a solution 
of the inverse problem follows from Theorem 15.31 It follows from Theorem 
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also, that the matrices {sj-k\kj=o generated by the Weyl functions are 
positive definite. 

Hence, it remaines to show that all the functions such that (I6.2ip holds 
and {sj-k}kj=o > (r > 0) are Weyl functions. Indeed, fixing such a matrix 
function ip, we get a sequence of matrices S(r) = {sj-jfe}& j=o > ®- Therefore 
we get a sequence of the transfer matrix functions wa{t, A) of the form ( 15.41(1 . 
where 



n(r) 



a 

a + s-i 



+ s. 



(6.22) 



and (15.241) holds. Taking into account formulas (15.411) . (15.47(1 and ( 15.43(1 . 
( 15.46(1 . we obtain matrix functions 



/?(r)*/3(r) :=U(r)*S{r)- 1 P*(P3(r)- 1 P*) PS^y 1 ^) (r > 0). (6.23) 
In view of the matrix identity (15.241) we have 



PS , (r) -1 n(r)Jn(r)*S , (r)- 1 .P* = -iP(s(r)- 1 A(r) - A^yS^r)- 1 ^ 



P5(r)- x P*. 



P* 



(6.24) 



In other words /3(r) satisfies the second relation in (I5.13p . Therefore formulas 
C r = 2Kf5(r)*/3(r)K — j define a system of our class on the semiaxis. In fact, 
this construction coincides with (15 . 19(1 and the function ip is the Weyl function 
of our system. Indeed, similar to the proof of Theorem 15.31 we derive from 
(I6.23P the equality (I5.4U() . Compare now Definition 15.41 and Theorem 16.21 to 
see that ip G Af(N) for any N. According to Theorems 16.41 and 16.51 it means 
that ip is the Weyl function. ■ 

Finally, consider the upper halfplane and define holomorphic Weyl func- 
tions in C + via relations (14.26(1 and (14.27(1 too. Put 



A7(iV) : = {ip(X)* : ip e Af(N)}. 



(6.25) 
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Remark 6.7 Similar to the proof that p = p, where tp and tp are given by 
( Ifi.ffj) and 116.6}) . respectively, one can show that the set J7(N) consists of lin- 
ear fractional transformations A5.10\) . where the pairs R, Q are meromorphic 
in C +; are well defined at A = i, and have the property 

R(X)*R(X) + Q(X)*Q(X) > 0, R(X)*R(X) > Q(X)*Q(X), X e C+. 

(6.26) 

In view of Remark 16.71 we obtain in C + the analog of Theorem 16.51 and the 
proof is similar. 

Theorem 6.8 Let system U.fy) be given on the semiaxis k > and let ma- 
trices C k satisfy A5.13\) . Then, the matrix function <£>oo(A)* = ("1^=0 -^(^) ^ s 
the unique Weyl function in C + of system on the semiaxis. 

Proof. Substitute <p 00(A)* instead of p 00(A) into (I6.16P and take into account 
( I6.26P to derive 



[^oo(A) I P ][J — q(X) r+1 KW r+1 (Xy 3 W r+1 (X)K* 



-Z^oo(A)* 



<2^oo(A)-^oo(A)*). (6.27) 

Now, inequalities (16.181) and (I4.26P are straightforward, i.e., ^^(A)* is a Weyl 
function. 

Instead of (16.191) we use the inequality 

q(X) k W k (XyC k W k (X) > -q(X) k W k (Xy 3 W k (X) > -j (X e C+), (6.28) 



which yields inequality 



J2[i P - i P ]qW k Kw k (xyc k w k (x)K* 



k=0 



-Ip 



oo. (6.29) 



The uniqueness of the Weyl function follows from (16.291) 



Theorem 16.81 for the scalar case p = 1 has been proved earlier in [TSJ, [23] (see 
also Theorem 3.2.11 [38J). 
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